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A Rayleigh-Schrodinger perturbation theory based on the Gaussian wavefunctional is constructed. 
The method can be used for calculating the energies of both the vacuum and the excited states. A 
model calculation is carried out for the vacuum state of the X4> 4 field theory. 
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(N ' I. INTRODUCTION 

c '■ 

1 & \ The Gaussian wave-functional approach (GWFA) [Q, or Gaussian approximation, has become a powerful and 

important tool to extract non-perturbative result of quantum field theory , finite temperature field theory |5| , and 
^ condensed matter systems [||J7] since Stevenson's advocation of the method about two decades ago ||. The Gaussian 
■effective potential (GEP) obtained from the GWFA provides a good starting point for further investigations of various 
I 'systems p-|l7|. Moreover, this approximation can also be used for realizing some novel ideas However, the GWFA 
.is essentially a variational approximation, and hence improvement of the obtained result is not straightorward. So 
^T) 'far, there exist mainly two methods to improve the GWFA. One way is to continue using the variational method with 
<^> more elaborate, non-Gaussian trial wavefunctionals. For example, Kummel and his collabrators developed the coupled 
~~ * 'duster method and obtained results beyond the Gaussian approximation of the \<j) 4 and <\P models ||; Ritschel and 
his collaborators constructed a non-Gaussian trial wave functionals through a nonlinear canonical transformation 
JToj] ; In order to investigate (3 -I- l)-dimentional Xdfi field theory, Yotsuyanagi proposed an improved scheme of the 
GWFA by adopting a BCS-type wavefunctional frLl|j . Another way of improveing the GWFA is to use appropriate 
expansions which give the GEP in their lowest order. In this aspect, Okopihska developed an optimized expansion 
'method to calculate the generating functional with the Euclidean formalism ; In the same Euclidean formalism, 
■ .Stancu and Stevenson formulated a slightly different expansion scheme and calculated the post-GEP in the spirit of 
"^■the background-field method ; Based on the GEP, Cea proposed a generalized GEP with a variational basis and 



o 



X 



carried out the calculation with the help of the standard perturbation technique in quantum field theory |l4| ] ; In the 
late 1990s, within the Minkowski formalism, one of the authors (Yee) and his collabrators developed the background 
field method to give an expansion of the effective action around the Gaussian approximate results |l5| ]. Recently, 
in order to calculate the partition function of a fermionic system, two of the authors (Kim and Nahm) and their 
collaborator proposed a variational perturbation scheme based on the functional integral without resorting to the 
background field method Additionally, Solovtsov et al. proposed a kind of variational perturbation theory to 
calculate the effective potential jl?]]. All these investigations improve the GEP or Gaussian-approximate result with 
miscellaneous degrees of success. Nevertheless, further investigations are needed still to achieve better approximations. 

We note that Rayleigh-Schrodinger perturbation theory is one of the basic techniques in nonrelativistic quantum 
mechanics [[l9] . It is also widely used in statistical mechanics |2Cj and, recently, generalized to quantum field theory 
pl| . Conventionally, Rayleigh-Schrodinger perturbation theory is based on an exactly solvable part of the Hamiltonian 
of a system. Recently, two of the authors (Kim and Nahm) and their collabarators calculated the energies of an 
anharmonic oscillator by combining Rayleigh-Schrodinger perturbation theory and variational method and obtained 
satisfactory results p2[ . In this paper, we develop a variational Rayleigh-Schrodinger perturbation theory based on 
the Gaussian wave-functional approach (RSPTGA) in quantum field theory. Different from the schemes mentioned in 
the last paragraph, RSPTGA can be used to calculate the energies of both the vacuum and excited states. Applying 
RSPTGA to A0 4 model, we calculate the vacuum state energy up to the third order. The result improves the GEP 
substantially and indicates a fast convergence. Comparing with the result in Ref. |Q, our result up to second order 
is shown to introduce an additional term. 
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Next section, we construct a quasi-free-field eigenstate set based on the GWFA, and develop RSPTGA based on 
the eigenset. In Sect. Ill, the vacuum state energy of the \(j> field theory will be calculated using RSPTGA. Finally, 
physical implication of the present results will be discussed. 



II. RAYLEIGH-SCHRODINGER PERTURBATION THEORY BASED ON THE GWFA 

In this section, we briefly introduce the GWFA, and construct a quasi-free-field eigenstate set. Based on this 
eigenstate set, we construct RSPTGA. 

We consider a model with the Lagrangian density 



C = -d^ x d^ x - V{4> x ) , (1) 



where x = (x 1 , x 2 , ■ ■ ■ , x D ) represents a position in D-dimcnsional space, (f> x = <j>(x) is the field at x, and the potential 
V(<p x ) has a Fourier reprensentation in a sense of tempered distributions p3fl . Many model potentials, such as 
polynomial models, sine-Gordon and sinh-Gordon models, possess this property. According to Ref. (J] (JPA), the 
GWFA produces the best trial vacuum wavefunctional 



(o) 1 f 

|0 > = A/"exp{-- / (<f) x - <p)f x y{<f>y ~ if)} 



(2) 



where, Af is the normalization constant (i.e., < 1 >= 1), f xy = J d D xd D y, and f xy — J d D pf(p)e lp< - x ^ with 
p = (p^p 2 , • • • ,p D ). The classical constant <p is equal to the Gaussian- vacuum expectation value of <fi x , 



(o) (o) 

<p=< 010,10 > . (3) 



Here, 



f(p) = vV+M 2 (^) , (4) 

and 



da 



H 2 &)= I ^e-^y (2) (-VMR^)) + ^), (5) 



oo 



where V^(z) = = d^(i q )ny^ e iqz (y(q) is the Fourier representation of V(zj) and I n {Q) 

I fffr ( ^+q?' ■ Thc GEP of thc s y stcm ' Eq-W' is s iven h y 

1 . . T . , f°° da 



V G (p) = -/o(A*)-^/i(M)+ / -^e-^V(-^hUT)+ V ) (6) 



2 4 J_ 00 2 y /^ V 2 

with fi chosen from the three possibilities: solution of Eq.(5), fj, = and /i — > oo @(JPG) (Hereafter, we will write 
(j.(<p) as fj, for simplicity except for special cases). We note that when Vg(</ ? ) has the absolute minimum at ipo, fJ-ifo) 
becomes just the physical mass and Vg(</?o) represents the vacuum state energy. The symmetry of the vacuum can 
be discussed using Eqs.(3)— (6). 

For the Gaussian vacuum, Eq.(2), one can construct the following annihilation and creation operators [Q (ZPC): 

A f (p) = (^7^) 1/2 I e^lfm* - 9) + at] (7) 

and 

4 (P) - ( 2(2 ^ /(p) ) 1/2 I e ipX lf(P)(4>* -V)- 01.] (8) 

with the relations [Af(p), A ^{p')] — 5{p' — p) and Af(p)\Q >= 0. Here, II X = —ig^- is the canonical conjugate 
operator to <f> x with the commutation relation, [<fi x , IL y ] = iS(x — y). Based on these operators, one can construct the 
quasi-free-field Hamiltonian, 
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H = J dpf(p)A\( P )A f (p) = |[in* + \{d x 4> x f + \f{<t> x - <p) 2 - ll (ji)] . (9) 

The Gaussian vacuum, Eq.(2), is the ground state of Hq with zero energy eigenvalue Eq°\ The excited states of Hq 
are [fHJ 

;o) i ™ . (o) 

>=— JJ4(pj)|0 >, n=l,2,---,oo (10) 



(0) 

In 

with the corresponding energy eigenvalue 



i4 0) = £/(?*)• (ii) 



Note that Eq.(2) is not a naive vacuum, since it contains information on the interacting system, Eq.(l). Obviously, 

(0) (0) , (0) (0) 

the wavefunctionals \n > and |0 > are orthogonal and normalized to < n|n >— X)p(n) II fc=i $(Pk ~ Pik) ( nere 

(o) 

Pi(n) represents a permutation of the set {ik} = {1, 2, • • • , n} and the summation is over all Pj(n)s). |n > describes 

(o) (o) 

a n-particle state with the continuous momenta pi,P2, • ■ • ,Pn- |0 > and |n > with n = 1,2, • • • ,oo constitute the 
complete set for Hq, which we call quasi- free-field complete set. 

Based on the above complete set, one can readily apply the conventional Rayleigh-Schodinger perturbation technique 
to calculate the energy of the system, Eq.(l). In order to do so, we write the Hamiltonian of the system Eq.(l) as 
H = H + H T = H + (H - H ) with 



j[-\f{<t>*-v? + \h{») + V(<i> x )) . 



Hi = / h^Ofe - + r/o(/i) + V{<f> x )] . (12) 



Following the Rayleigh-Schodinger perturbation procedure [ [19pl] ], one can obtain the wavefunctionals and energies 
for the vacuum and excited states respectively, 

\n >= £[Q« ^(E n - 4 0) - Hi)] l \n > , (13) 

i=o H - E n 

and 

M<P) = E( n ] + E < *\Hi[Q» ~ — (oyC-^n " 4° 3 " A"/)]'!" > (14) 

J=0 #0 - E n 

with Q„ = E,y„ J d D Pirf D P2 • • • d 2 ^ |j >< i|. 

For the case of n = and Z = 0, Eq.(14) gives the vacuum enegy up to the first order Eq = E^ + Eq 1 ^ — J Vg(^>) 
which is just the product of the GEP and the space volume. Thus, employing Eq.(14) with n — 0, one can get the 
effective potential for the system, Eq.(l), 

*.(•)- $g (-) 

<owgo> 

<0|0> ' v ' 

which takes the GEP, Eq.(6), as the first-order approximation. Eq.(16) implies that $ should replace tp in the 
calculation of the effective potential. It is evident that <p in Eq.(3) is the zeroth order approximation, of <F 
When Eq.(15) is truncated at nth order, one should also truncate Eq.(16) at the same order |l2|, 13]] . 

In the case of n ^ 0, Eq.(14) is the excited-state effective potential of the system [||. Additionally, Eq.(14) becomes 
the vacuum and excited-state energies for the symmetric phase of the system if ip = is chosen in the scheme. In the 
next section, we apply RSPTGA to the vacuum state of the A0 4 field theory. 



with 
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III. APPLICATION TO FIELD THEORY 



In this section, we consider the potential, V{<j> x ) = ^m 2 <j) x + ^4>%, which was widely studied in connection with the 
Gaussian approximation [T^ ] . For this system, Eq. (5) gives rise to 

^ = m 2 + l -\^ + - A \h{n), (17) 

and from Eq.(6), one readily obtains 

V G M = imV + + \h{p) - ^I(V) , (18) 

which is just Eq.(2.22) in Ref. Jl|l(PRD) Q To obtain the effective potential of the A0 4 field theory up to a given 
order according to Eq.(15), we first calculate the following matrix elements of Hi : 

< n|ffz|n > = < 0\Hj\0 >< n\n > E II A ^ "^J 

1 ^ {P ;j (n-2)} fe=l 



and 



with 



%L-i) -Pi»)[/(Pi (n _ 1) )/CP<(„-i))/(Pi»)] 2 

P,,(n-l)fc=l ' 



(20) 



(°) I rr I „ (0) 1 A V- TT A / / ^ 

< nlff/h-2 > = ; -pr > IT A(p h - p' 4 

■ (5 (^(„-2) _ P*r»-2 _ Pi(n-l) ~Pi„)[/(pi- (n _ a ))/(Pi»- a )/(Pi(„-l))/(Pi»)]~' , (21) 

^n!(n- 3)! 2 V2(2vr) c p ^ 3) fe=1 

•*(Pa ( „- a) +Pi ( „-i) + jO[/(Pi ( «- a ))/(Pi ( »-i))/(PiJ]~* (22) 



(o) ,„ (o) 1 A „ " 4 

•<%(„-3) +Pin-2 + Pi(„-i) + P»n ) [/ (P»(„- 3 ) ) /(Pi(„-a) )/(Pi(»-i) )/(Pi» )] ~ ' > (23) 



0, for fc < 

A(p iik -pJJ=^ 1, for fc = 

${Pt k -p' 3k ), for fc>0 



x our notation /„(» is different from Eq.(2.21) in Ref. @(PRD). 
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(0) 

Here, the index i& £ {1, 2, • • • , n} with fc = 1, 2, • • • , n corresponds to \n > and jk S {1, 2, ■ • • , n } with fe = 1, 2, • • • , n 

(o) 

to < n \. Pij (I) represents a given permutation of I momenta ,Pi 2 , ■ ■ ■ , paired respectively with , , • • • , p- , 
and X)p (;) ^ s over all different Pij(l)s. For any Py(Z), ii,i2> ■ " »*i are different from one another, and so are 
ji , j2 , ■ ■ ■ Employing the above matrix elements, a straightforward, yet lengthy calculation according to Eq.(15) 
gives the effective potential of the A</> 4 field theory up to the third order as 

Tpiii 

11 9/2 ^ 2\2 ^ A 2 2 B \ 2 

A + A! + A 2 A 2 2 2 A 2 2_B 1 _f5 2 _A 3 2 C A 3 
+ 48 ^ ^ " 3^ } + + 512^ (24) 



with 



A = I 4> (ag^ [/i(»)/i(y)/i(» + y)]- 1 !/!^) + AG/) + A(* + y)]" 1 , 
^ = / (^^[/iW/iW/i^ + y)rMi + AW + A (2/) + AO* + vT 1 , 
A 2 = I + yrUM + A(y) + A(* + y)]- 1 

■ii+/iW + /i(!/) + /i(x+»)r i , 

s = / (2 ^ l/iWAM/iW/i^ + y + 

•[/iW + AG/) + A(*) + + y + z)]" 1 , 

fl i = / ^(^^I/iW/iW/iW/i^ + v)h(* + y + Z T' 

■ifi(x) + AG/) +A(x + yT'ifiix) + h{y) + h{z) + Mx + y + z)}' 1 

(2^ (2% (2^ [h{x)h{y)h{z)h{x + y)h(x + z)]- 1 
■[fi(x) + h(y) + h(x + yT^hix) + fi(z) + h(x + z)]- 1 
C = I 4> ($» { ^ \h{x)h{y)h{z)h{u;)h{x + y + z)h{x + y + c)]- 1 

• [A (x) + fi (y) +fi(z) + fi(x + y + z)} - 1 [h (x) + ft(y) + AM + fi(x + y + w)Y 



and fi(w) = vi + w 2 . Here, we take ip as the zeroth order approximation of In Eq.(24), the second line is the 
second order correction and the third line the third order correction to the GEP. 

From Eq.(24), one can see that after Vg(v) is renormalized, its corrections will be finite and, accordingly, a further 
renormalization procedure is not needed for higer order corrections. In the case of (0+1) dimensions, our second-order 
result with tp = is consistent with Eq.(14) in Ref. [ p2| . Furthermore, numerical calculation for the case of (1+1) 
dimensions indicates that the effective potential with the second-order correction predicts existence of a second-order 
phase transition. We also note that the second-order correction improves GEP substantially, and the third order 
correction is vanishingly small. This can be seen by comparing the coeficients of the first, second and third terms in 
the second line of Eq.(24) with corresponding ones in the third line, respectively. This indicates that RSPTGA has a 
fast convergence in general. 



IV. DISCUSSION AND CONCLUSION 



In this paper, a Rayleigh-Schrodinger perturbation theory based on the GWFA within the framework of quantum 
field theory is proposed. Since the theory is based on the Gaussian approximation, it provides a systematical tool for 
controlling the Gaussian approximation. It can be used not only for calculating the effective potential but also for 
considering excited states. When one is interested in symmetric phase, the vacuum and excited-state energies can be 
calculated beyond the Gaussian approximation by RSPTGA using ip = 0. Application of RSPTGA to the A</> 4 field 
theory shows that it can improve the GEP substantially with a fast convergence. 

We note that RSPTGA predicts existence of a second-order phase transition in the (l+l)-dimensional Xtfi 4 field 
theory, although the critical coupling is very small. It may be attributed to the fact that for the second-order case 
we approximated $ =< 1 1 > L_ 2 as tp. We also note that the second-order result of Eq.(24) has an additional 
term — ^-i^ip 2 ^ 2 — ^<p 2 ) 2 which does not appear in Ref. (li). In fact, Eq.(12) gives, for the A0 4 field theory, 

E I = / {V G (p) + ^ 2 - ^ 2 ) : {4> x - p) : +±tp : (<fe - <^) 3 : : (<k - ■}, (25) 
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which is equal to Eq.(4.19) in Ref. jl4| (PRD) except the J Vc{ip) term, where the colons mean normal ordering with 
respect to the Gaussian vacuum, Eq.(2). The additional term in Eq.(24) arises from the linear term </?(/Lt 2 — ft/? 2 ) : 
(<f) x — if) :, and disappears when the constraint Eq.(5.18) in Ref. 113] (PRD) is adopted in our scheme; that is, when 
$ in Eq.(16) is taken as (p. However, we note that tp is simply the Gaussian- vacuum average value of the field 
operator (see Eq.(3)) and $ is not equal to ip when higher order contributions are calculated Jl(],|l^,|l3|] . Moreover, 
the variational procedure which led to the GEP, Eq.(18) produces the extremum condition, Eq.(17) (i.e., Eq.(2.19) 
in Ref. jl4j (PRD)), and Eq.(5.18) in Ref. [H (PRD) at any truncated order is not compatible with this extremum 
condition. 

In closing the paper, we like to point out that it is straightforward to generalize RSPTGA to other cases, such 
as a Fermion field system. In fact, a Rayleigh-Schrodinger perturbation technique based on the variational results 
has been applied to a polaron problem [£3|. Moreover, one of the authors (Lu) developed the GWFA in thermofield 
dynamics Q ||(Lu). Based on it, it is possible to develop RSPTGA within the framework of thermofield dynamics 
which will be useful for treating finite temperature cases. Finally, instead of performing the variational procedure in 
the GWFA as described above, the extremization process with respect to /i can be carried out after truncating the 
series of Eq.(14) at some given order jL2|. This procedure will lead to a slightly different variaion of RSPTGA ||. It 
may have its own advantages or peculiarities over the scheme developed here. 
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